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Abstract 

The macroscopic quantum tunneling of a planar domain wall in a ferromag- 
netic metal is studied based on the Hubbard model. It is found that the ohmic 
dissipation is present even at zero temperature due to the gapless Stoner exci- 
tation, which is the crucial difference from the case of the insulating magnet. 
The dissipative effect is calculated as a function of width of the wall and is 
shown to be effective in a thin wall and in a weak ferromagnet. The results 
are discussed in the light of recent experiments on ferromagnets with strong 
anisotropy. 
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In recent years, owing mainly to the development of technology in mesoscopic physics, 
there has been growing interest in macroscopic quantum tunneling (MQT) in magnetic 
systems e.g., the magnetization reversal in small grains 0, the quantum nucleation of 
a domain p|, and the quantum depinning of a domain wall via MQT These studies 
are mainly in ferromagnets, but recently a magnetization reversal due to MQT has been 
observed in antiferromagnetic particles of horse spleen ferritin [^. In the case of the quantum 
depinning of a domain wall pinned by defects, the position of the wall at the pinning center 
becomes metastable in the external magnetic field, and if the barrier height is low enough, 
the position tunnels out of the local minimum. This problem was studied theoretically by 
Stamp ^ for the case of an insulating magnet. The tunneling rate was expressed in terms 
of macroscopic variables, and was shown to be large enough to be observed even for a large 
wall with about 10^° spins. As sources of dissipation, which is shown to be important by 
the seminal paper by Caldeira and Leggett P, Stamp considered magnons and phonons, 
but the effects turn out to be negligible, since magnon has a gap and coupling to phonon is 
weak. Consequently it has been concluded that the tunneling rate is not essentially affected 
by dissipation in insulators. 

Experiments on MQT in magnetic system, however, have been carried out in metallic 
ferromagnets. In metals, in contrast to the case of insulators, there is a gapless excitation 
of spin flip, and hence dissipation from conduction electrons must be very important. Con- 
sequently the quantum motion of the wall in metals should be quite different from that in 
insulators. In this paper, we will investigate theoretically the dissipative effect on MQT of 
a domain wall in a ferromagnetic metal based on an itinerant electron model. An important 
and interesting feature of the itinerant system is that the electron, which supports magne- 
tization, works also as a source of dissipation in the dynamical motion of the magnetization 
itself. Our analysis is based on the Hubbard model in the continuum. The calculation 
is carried out at zero temperature, since we are interested only in the quantum tunneling 
present at low temperature. 

The Lagrangian in the imagnary time path integral is given by 
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(1) 




where Cxo- is an electron operator at site x with spin a{= ±), rixo- = c\.„c^fj and U is 
the Coulomb repulsion. The band energy is Ck = k^/ (2m) — ep with the fermi energy ep. 
The Coulomb repulsion term will be rewritten by introducing a Hubbard- Stratonovich field 
representing the magnetization; Mx = MxHx, where Mx =< (c^<tc)x > Hx with Hx being a 
slowly varying unit vector which describes the direction of magnetization. The magnitude 
of magnetization is assumed as space-time independent, Mx = M. Hence only nx remains 
as the relevant degree of freedom. 

The spatial variation of nx accompanied with a domain wall is assumed to be much 
slower compared to the inverse fermi momentum of the electron k^^. For the analysis of 
such a slowly varying field, a locally rotated frame of electron is convenient such that the 
z-axis of the electron is chosen in the direction of the local magnetization vector nx. The 
electron operator Ox^ in the new frame is related to the original Cxa as 



where the polar coordinates (^^x(t), 0x(t)) parametrize the direction of nx(r). The electron 
Oxer is polarized uniformly with the energy Cko- = k^/2m — aUM — ep. As a price of this 
transformation, there arises from the kinetic term c^c+ |Vcp/(2m) an additional term ifjnt 
that describes the interaction of electrons with spatial variation of the magnetization vector 
P). This interaction H[^t is small and of the order of 0(A;f|A)^^, where A is the domain 
wall thickness, /cpt is the fermi momentum of the majority spin, and hence can be treated 
perturbatively. Our following results are valid for Afcpt ^ 1- 

The integration over the electron degrees of freedom leads to the effective action for the 
magnetization as Ses = — trln(9^ + ek<j) + PJ2x{U/2)]['P + AS. The first two terms are the 
mean field action for a ferromagnet which determines the magnetization M. The dynamics 
of {9, (p) is described by AS*, which is expressed in terms of correlation functions of electron. 
This term is decomposed into two parts, that is local and non-local in time, respectively, as 
AS* = AS'ioc + A^dis- The local part AS'ioc determines the dynamics of magnetization vector. 



a- 



acos(e/2)cx<x + e"''"'' sin(^/2)cx,-<. 



(2) 
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and the non-local part AS'dis represents the dissipative effect due to conduction electrons on 
the motion of the magnetization vector. 

Up to the lowest order in dr and V, the local part ASioc turns out to be formally the 
same as the ferromagnetic Heisenberg model with spin 5* = M/2 whose Lagrangian is 
given by 

S JS^ 



/s / s 

rf^x 2-0(1 - cos^) + (iyef + sin^ ^(V0) 



(3) 



The exchange coupling or the spin wave stiffness constant is expressed by the parameters of 



where 



the original Hubbard model as J = {n/marM^) 1 — ((/cpt ~ kp^)a^) / {?>QT[^mnU M) 
n is the electron number per site, /cpo- = (2m(eF + crUM)Y^'^ is the fermi momentum and 
a is the lattice constant. Hence, in the absence of the non-local term AS'dis, there is no 
formal difference between metallic and insulating ferromagnets, and the tunneling rate of 
the domain wall is determined on the same footing [Q. 

In order to incorperate the domain wall, the anisotropy energy with yz easy plane is 



introduced 10 



/fani = / f/'x (^-^S' cos' e + sin^ e cos^ 0) . (4) 

The Lagrangian Lh + -ffani has a planar domain wall centered at x = Q{t) and moving 
slowly as a classical solution; cos 6'(x,r) = tanh(x — Q{t)/X) and cos0(x, r) ~ iQ/c <^ 1 



with c = K±XSa^ where X = J/K is the width of the wall. This configuration is depicted 
in Fig. 

For the magnetic field H close to the coercive field He, i.e., {He — H) / H^ = e <^ 1, the 
potential for the wall coordinate Q is given by V{Q) = (l/2)M^a;Q(5^[l — [Q/QoY] where 
Myj = 2N / {Kx^X'a?) is the domain wall mass being the number of the spins in the wall. For 
this case of small e, the attempt frequency around the minimum is ujq ~ {fio{h'y)'' /a^)y/h^e^ 



and the width of the barrier is given by Qq = y3/2y/eX where he = Hd {fi'yS / a?) is the 
ratio of the coercive field to the magnetic moment per unit volume (/io is the magnetic 
peameability of free space and 7 is the gyromagnetic ratio). The actual value of attempt 



frequency is cuq — 5 x y^e^ (K) for the choice of a = 3 A, and in the present case, this 
is roughly the same as the crossover temperature T^o from the thermal activation to the 
quantum tunneling. The classical solution (bounce) of Q in the metastable potential V{Q) 
is given by Q{t) = Qo/cosh'^{uJoT/2), and the tunneling rate out of the local minimum is 
estimated by use of this bounce solution. For the case of the wall with the cross sectional 
area Na^/ A as shown in Fig. the rate Fq without dissipation is reduced to Fq = A exp{—B) 
where A ~ {fio{h'y)'^ /a^)N^/'^hl/'^e'^^^ =i W^N^/'^hl/^e'^/^ (Hz) and the exponent,^, is given 
as -B ^ Nh^J'^e'/'^. 0] Since B is proportional to Ne'/^ a small value of e is needed to 
observe the tunneling. 

Let us now look into the non-local action AS'dis, where the characteristic feature of the 
itinerant electron system is to be seen. For the case of a weak dissipation, this contribution is 
evaluated by use of the configuration of a domain wall obtained in the absence of dissipation. 
Up to V^, AS'dis is obtained as 

x|9<,(r) - 9<,(r')|" < J;(q)Jl(-q) > |.„„ (5) 

where Ja(q) (« = ±,2;) are the Fourier transform of the spin currents of the electron; 

= —i[{a)aaVa) — (VaV^a)] with a± = ax^ iCy, and 6'q = I]x ^~*''''^x- The dissipation 
does not result from the z-component Jz in the present case of a domain wall motion with 
V0 = 0. The expectation value of electron spin current < J+ J_ > in AS'dis is evaluated by 
the random phase approximation (RPA) in the background of uniform magnetization. |I2| 

After the analytic continuation to real frequency, ASdis is expressed by the imaginary 
part of the retarded correlation function < J+ J_ > \uj+io as |13| 



-/m2 



X r ^e--l---'llm < J^(q) J^(-q) > U+.o- (6) 

Jo TT 



The imaginary part is expanded in terms of w/ep as 
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Im < J^(q)J^(-q) > 



UJ 



7rlg|3 



+ Oiyfi^ /cpt — ^Fi < < ^Ft + ^Fi 



(7) 







otherwise. 



The term hnear in uj gives rise to the ohmic dissipation. It is seen from the restriction on 
q that the ohmic dissipation arises from the Stoner excitation, which is a gapless excitation 
of spin flip across the fermi energy. 

By the expression of the domain wall configuration, the non-local part of the effective 
action is reduced to 



electrons and the wall, and because of this factor, the momentum integration is domi- 
nated by g < A^^. The time integral is estimated by approximating the bounce solution as 
Q(t) ~ Qi^Q{uj^^ — |r|) and by introducing a short time cutoff of uJq^ for the relative time 
(r — r') Noting gQo q\\ft ^ 1, the sine function in Eq. (|) can be replaced by its 
argument and the action is evaluated to be AS'dis = TiNe where the factor, e, is due to the 
smallness of the squared tunnel distance Q\. Here the strength of dissipation, r], is 



For a thick wall A(A;ft "~ ^fj,) ^ 1, oc exp[— 7rA(A;Ft "~ ^Fi)] and then the dissipation 
is negligible. On the other hand, ri can be large if (/cpt ~ ^Fi)A < 1. This condition is 
compatible with that of slow spatial variation A/cpt ^ 1 for a wall with moderate thickness 
in a weak ferromagnet and for a thin wall in a stronger ferromagnet. The strength ri is plotted 
as a function of (A/a) in Fig. for three different values of 5 = (A;f| — A;Fi)/(fcFT + ^fj,) 
with (/cpt + ^Fi)a. = 6 which may represent the case of an iron. The dissipation is larger 
for weaker magnet (smaller S). (For a complete ferromagnet, fep^ vanishes and the ohmic 
dissipation disappears.) It is seen that r/ can be of the order 0.1 for a wall with thickness 





The form factor of the wall, 1/ cosh^(7rAg/2), represents the effective coupling between 
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a few times the lattice spacing with 5 < 0.1. In the presence of dissipation, the tunnehng 
rate is reduced to be F = Aexp[—{B + A^dis)] = Fq exp(— r^A^e). Because of the different 
e-dependence of B and ASdis, the ratio A^disZ-B = rjh'^^^'^e^^^^ is much larger than unity for 
the case of small e we are interested in, and in particular for a thin wall [h^ is usually small, 
e.g.,~ 10~^). Consequently the tunneling rate is predominantly determined by dissipation 
in such cases. The tunneling rate F is shown in Fig. |] for the case of insulator [r] = 0) 
and the typical case of a metal {t] = 0.1) by the broken and solid lines, respectively for a 
choice of he = 10~^. In this figure, the number of spins is taken either = 10^ or 10^. 
The value N = 10^ corresponds, for instance, to a wall with thickness of about lOA and the 
area of 200Ax200A. The tunneling rate is seen to be much smaller in metals than that in 
insulators. 

We have neglected the effect of magnetic field on electronic states. This is justified as long 
as UM ^ 'yH. In experimental situations with the magnetic field of < IT and U ~ lOeV, 
this condition reduces to M > 10~^ in unit of the Bohr magneton, which is easy to satisfy. 
However, in order to discuss the case of very small M, the fluctuation of the magnitude Mx 
around the mean field value must also be included. 

The contributions of higher order in H^^t are smaller than that of the second order we have 
calculated; for the potential renormalization by the order of (kpi^X)'"^ and for the dissipative 
effect by (/cf|A)~^ or e. 

In Eq. d^) we have taken account of only the ohmic dissipation. The super-ohmic 
contributions, which are of higher orders of {uj/ep) in Eq. (0), are smaller than the ohmic 
one by a factor of (cjo/ep)^ <^ 1 and hence are negligible. On the other hand, a contribution 
from the magnon pole, which has not been taken into account in the correlation function 
< J J >, is calculated from 

Im < J^(q)J-_(-q) > |SViJ ^ nM'{Jmqf6{u - u,) (10) 

where Ug = Aq + is the magnon energy with the anisotropy gap Aq. This 

pole leads to super-ohmic dissipation, whose strength, 

^(pole)^ is evaluated as r^'P" = 
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(28/5)M(Ao/ct;o) exp(— Ao/c<Jo)- Since experiments are usually carried out in highly 
anisotropic materials with Aq/cjo — 10, this contribution is very small compared to the 
ohmic dissipation for the case of a thin wall. 

The present metallic case, where the ohmic dissipation is present even at absolute zero, 
are in contrast with the insulating case. At finite temperatures, however, there are ohmic 
dissipations even in the latter case. Stamp calculated such ohmic dissipations from two- and 
three-magnon processes and found r]^^^^^ =(3/27rS')(l//?Ao)exp(— /3Ao). These processes 
corresponds to higher order contribution of i^int in our calculation. In contrast to the case 
of metals, the value of r/^™'*^) vanishes at T = and is very small at /3Ao ^ 1, hence the 
ratio of AS'dis/-B is negligibly small in insulators. 

In metals, eddy currents may influence MQT. An electric field is induced by Faraday's 
law from the change of magnetizations accompanied with the motion of the wall. This 
field produces the electric current and thus leads to the Joule heat of P = (/io^7/a^)^(5^/p 
per unit volume where p is the resistivity. By use of the specific heat C and the system 
size of L, the temperature rise due to the eddy current is expressed as AT = PL/{CQ). 
For p ~ 10^^[fim], C ~ lOfJ/Km^], and Q ~ l[m/s] with L = lOOA, it is estimated as 
AT ~ lOmK. This value is not negligible but effects associated with this temperature rise 
may be separated from the intrinsic effects in careful experiments. 

Our result shows a distinct difference between MQT of thin walls in metallic and insu- 
lating magnets. Unfortunately the experiments carried out so far appears not yet be able to 
observe dissipation due to itinerant electrons. In the experiment of a domain wall motion in 



a small particle of Tbo.5Ceo.5Fe2, MQT was observed below T^o — 0.6K [0. However, the 
width of the domain wall is about 30A and according to our result, 77 oc exp[— 7rA(/i;FT — ^Fi)], 
the dissipation from electron is negligible for such a thick wall. This may be the reason why 
the result of the crossover temperature Tco ~ 0.6K is roughly in agreement with the theory 
Q without dissipation. On the other hand, the domain wall in SmCo5 is very thin A ~ 12A, 
and our result suggests strong effect of dissipation, which will be interesting to observe. 
Experiments on bulk crystal of SmCo3.5Cui.5 with very thin walls (a few times a) have been 
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performed |T6|, although quantitative argument is not easy since many walls participate in 
these experiments. Even in the case of thick walls, the dissipative effect becomes large in 
weak ferromagnets, where the experiments, however, will not be easy because of small value 
of saturation magnetization M. 

MQT in disordered magnets has a new possiblity of observing a significant effect of 
sub-ohmic dissipation. In fact, as disorder is increased in a metallic magnet, the Anderson 
transition into an insulator will occur, and it was shown recently that near the transition 
the dissipation due to the conduction electron becomes sub-ohmic [|l3l . Disordered magnets 



may also be suitable for study of MQT because the eddy current becomes less important 
for larger resistivity. 

Our calculation is valid in the s-d model as well, where localized magnetic moment is 
due to d electron and the current is carried by s electron. 

In conclusion, we have studied the macroscopic quantum tunneling of a domain wall 
in a metallic ferromagnet on the basis of Hubbard model. The crucial difference from the 
case of an insulator is the presence of ohmic dissipation even at zero temperature due to 
the gapless Stoner excitation. The coupling of domain wall to electrons is effective only for 
momentum transfer of |g| < A^^, while Stoner excitation is gapless at the restricted region 
/cpt ~ k-pi < \(l\{< ^Ft + ^Fi)- Hence the effect is negligible for a thick domain wall in which 
experiments so far have been carried out. On the other hand, important effects of the ohmic 
dissipation are expected in thinner domain walls and in weak ferromagnets, which will be 
within the present experimental attainability. 
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FIGURES 



FIG. 1. Configuration of a planar domain wall. 
FIG. 2. Strength of ohmic dissipation rj given by Eq. (^) as a function of the width of the 

wall A/a, a being a lattice constant, for three choices of 6 = (/cpt ~ ^Fj,)/(^Ft + = 0.05,0.1 

and 0.2 with (kp^ + kpi)a = 6.0. 

FIG. 3. The tunneling rate F for the insulating (rj = 0) (dashed line) and the metallic {rj = 0.1) 

(solid line) magnet as a function of e, He being the coercive field. Number of spin is = 10^ and 

108. 
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